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After a review of temperature-dependent experimental x-ray absorption near-edge structure 
(XANES) and related theoretical developments, we present the Al ii'-edge XANES spectra of corun- 
dum and beryl for temperature ranging from 300K to 930K. These experimental results provide a 
first evidence of the role of thermal fluctuation in XANES at the Al A'-edge especially in the pre- 
edge region. The study is carried out by polarized XANES measurements of single crystals. For 
any orientation of the sample with respect to the x-ray beam, the pre-edge peak grows and shifts 
to lower energy with temperature. In addition temperature induces modifications in the position 
and intensities of the main XANES features. First-principles DFT calculations are performed for 
both compounds. They show that the pre-edge peak originates from forbidden ls^-3s transitions 
induced by vibrations. Three existing theoretical models are used to take vibrations into account in 
the absorption cross section calculations: i) an average of the XANES spectra over the thermal dis- 
placements of the absorbing atom around its equilibrium position, ii) a method based on the crude 
Born-Oppenheimer approximation where only the initial state is averaged over thermal displace- 
ments, iii) a convolution of the spectra obtained for the atoms at the equilibrium positions with an 
approximate phonon spectral function. The theoretical spectra so obtained permit to qualitatively 
understand the origin of the spectral modifications induced by temperature. However the correct 
treatment of thermal fluctuation in XANES spectroscopy requires more sophisticated theoretical 
tools. 
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I. INTRODUCTION AND STATE OF THE ART 



X ray absorption near-edge structure (XANES) spec- 
troscopy is a powerful technique to probe the empty 
states in solids and to determine the local structure 
around a selected atom. The interpretation of K- 
edge XANES spectra is not straightforward and often 
requires the use of simulation tools, which are tradi- 
tionally based on the density functional theory (DFT). 
XANES calculations for inorganic solids usually consider 
the atoms at fixed positions, even if in reality atoms 
are subjected to quantum thermal fluctuations that re- 
duce to the zero point motion at T=OK. Recently, it has 
been theoretically shown that vibrations could have a 
spectacular effect at the Al K-edge, allowing Is to 3s 
dipole transitions.^ Such forbidden transitions were pre- 
viously invoked to be responsible for the pre-edge feature 
that occurs at the Al i^-edge in several alumino-silicate 
minerals^ and at the Si K-edge in silicon diphosphatef2, 
Moreover, Ankudinov and Rehr stated that atomic dis- 
placements of order of typical Debye- Waller factors could 
reveal some forbidden transitions.^ Consequently, the 
fact that vibrations may be able to induce additional 
peaks (in the pre-edge or at higher energies) suggests 
that XANES could be considered as a relevant probe of 
quantum thermal fluctuations. It is already the case in 
organic molecules, the vibronic fine structure of which 
can be observed in XANES or NEXAFS (Near-Edge X- 
ray Absorption Fine Structure) spectra. For instance 
in the last ten years, XANES spectroscopy has revealed 



the vibronic structure of naphtalene^^ biphenyl^ and even 
larger molecules such as NTCDAi at the C K-edge, the 
vibronic structure of acetonitrile^ and acrylonitrile^ at 
the N K-edge, and more recently the one of halogenated 
acenaphthenequinoneiS at the O K-edge. The vibronic 
component of the lowest core-excited state of OH and 
OD was observed at the O i^-edge as wellJA 

In this paper, we focus on the thermal effects in two 
inorganic crystalline solids at the Al K-edge: corundum 
(a-AhOs) and beryl (BcaAlaSieOis). Polarized XANES 
spectra have been measured on single-crystals for tem- 
peratures ranging from 300K to 930K. Our purpose is to 
show that thermal fluctuations induce substantial spec- 
tral modifications, especially in the pre-edge region. Very 
few studies were carried out in inorganic solids towards 
the same goal. As far as we know, thermal effects mea- 
sured in XANES spectroscopy mainly concern the Ti K 
pre-edge region in oxides. Such a temperature depen- 
dence was first observed by Durmeyer et al^ in Ti02 
(rutile), Li4/3Ti5/304 and LiTi204. The main pre-edge 
peak was found to grow with temperature, and to slightly 
shift to lower energy for the first two compounds. Later, 
the temperature dependence of the polarized pre-edge 
structure at the Ti K-edge in rutile was measured at 
low and room temperature. L3, Thermal effects were essen- 
tially visible through an increase of intensity of the first 
two pre-edge peaks. The temperature dependence of the 
electric dipole and quadrupole contributions of the Ti K- 
edge in rutile was investigated in the 6K-698K range by 
Collins and Dmitrienko in 2010 at the Diamond Light 



Source. "'^'^ Nozawa et alJ^ measured the temperature de- 
pendence of the Ti JsT-edge XANES in SrTiOg from 15K 
to 300K, showing an increase of the first two pre-edge 
peaks along with a shght shift to lower energy of the sec- 
ond pre-edge peak and of the absorption edge. These 
effects were attributed to random thermal vibrations of 
Ti^ 

Apart from these last examples at the Ti K-edge, most 
of temperature-dependent XANES measurements were 
carried out in order to study phase transition mecan- 
isms or to determine the local structure below and above 
the transition temperature. For instance in some per- 
ovskite crystals, significant spectral changes in the K 
pre-edge region are induced by temperature, and are 
interpreted in relation to phase transition. The liter- 
ature reports temperature-dependent x-ray absorption 
spectra (including pre-edge) in perovskite titanates at 
the Ti K-edger^"— in perovskite manganites at the 
Mn K-edge^^^ and at the O K-edger^^ in zirconates 
at the Zr i^-edge^ niobiates at the Nb K-edge^^^^ 
in the Lai_a;Srj;Fe03_5 system at the Fe K-edge^^ 
in a La(Fe,Ni)03 solid oxide fuel cell cathode at the 
O K-edger^ and more recently in the Pro.sCao.sCoOs 
cobaltite at the Ca i2,3-edges and the Pr L^ and M4^5- 
edgesi^ Phase transition studies using temperature- 
dependent XANES are not restricted to the perovskite 
structure. In particular Cu i^-edge and La L-edges 
spectra in La2-a:Sra:Cu04 cuprates were recorded at 
room temperature and at 78Ki^ An angular dependence 
XANES study was performed at the V K-edge in VO2 for 
temperatures lower and higher than the metal-insulator 
transition temperature at 68° C;"^^ Several studies re- 
ported in literature deal with oxide glass and melt struc- 
tures. For instance, the local structural environment of 
Ti in Na-, K-, and Ca-titanosilicate glasses and melts was 
determined by Ti K-edge x-ray absorption spectroscopy 
at temperatures ranging from 293-1650 Ki^ The Zr K- 
edge in a Zr02-MgO-Al203-Si02 glass was measured for 
various temperatures around Tg (1085K)<^ The crys- 
talline and melt structures of AI2O3 and MgAl204 were 
investigated using temperature dependent XANES spec- 
troscopy at the Al and Mg isT-edgesi^ A similar study 
concerns the crystal, glass and melt structures in the 
CaO-MgO-Al203-Si02 system at the Al, Si and Ca K- 
edges4S Again at the Si i^-edge, the structure of Si02 
polymorphs were studied using XANES spectroscopy 
from room temperature to 2030K, above the melting 
point at 2000K4i Furthermore, at the oxygen ii'-edge in 
liquid water, variations were observed between the room 
temperature spectrum and the 90° C one and found to be 
related to a change of local environment of a significant 
amount of molecules J^ 

Temperature dependent x-ray absorption spectroscopy 
was also used to study chemical reactions. For instance, 
in catalysis, the local structure of Al in zeolites was ex- 
amined to understand the processes that take place dur- 
ing steam activation of these complex cagelike porous 
compounds. Hence Al .ftT-edge measurements in several 



zeolites were carried out in the range 300K-975K, in vac- 
uum or in a flow of helium saturated with water^^^— 
and in other workin g co nditions of the catalysis 41 In 
Refs. I43, [4^5 H, and [43 a pre-edge peak appears when 
temperature is increased. The presence of this Al K pre- 
edge peak at high temperature is interpreted as the sig- 
nature of three-fold coordinated Al in the zeolite struc- 
ture. Other temperature effects are observable at higher 
energies in the XANES region, and attributed to the 
transformation of a given amount of octahedrally coordi- 
nated Al to tetrahedrally coordinated A1.™4S, Recently, 
in a totally different field, the temperature dependence of 
Fe ii'-edge XANES spectra of FePt/Fe304 nanoparticles 
has been measured from 300K to 870K in order to un- 
derstand the annealing process, which enables the build- 
ing of magnetic nanocomposites with combined magnetic 
properties.— A strong temperature dependence was ob- 
served in Fe L2,3-edges spectra of Fe impurities in MgO 
thin films in the range 77K-500K; it was attributed to 
the thermal population of low- lying Fe 3d excited states, 
that are present due to the spin-orbit coupling 4^ At last, 
the valence state of Yb in YbC2 , investigated by Yb L3 
XANES spectroscopy at low and high temperatures, is 
found to be stable from 15K to 1123K, so providing new 
information about this long-known compoundi^ 

In all the references cited in the previous two para- 
graphs, the role of nuclear motion in the temperature- 
dependent XANES spectra was not investigated. A 
thorough study of thermal fluctuations could then bring 
new insights in the interpretation of the temperature- 
dependent spectral features whatever the purpose of the 
experiments is (catalysis, phase transition, etc.). 

The modeling of nuclear motion in the absorption 
cross section of solids is a challenging task. In partic- 
ular, the presence of phonon modes in crystalline solids 
generates a dynamic disorder, which has a completely 
different behavior from the static disorder that can be 
found in glasses. Far from the edge, in the extended x- 
ray absorption fine structure (EXAFS) region, vibrations 
are taken into account through a Debye- Waller factor 
exp (— 2fc^cr^)4ii^ Since this factor vanishes for /c ~ 0, it 
cannot reproduce well the thermal effects in the XANES 
region and especially in the pre-edge region. Fujikawa 
and co-workers showed in a series of papers of increasing 
sophisticationSi^ that the treatment of vibrations could 
be achieved by the convolution of the " phononless" x-ray 
absorption spectrum with the phonon spectral function. 
To our knowledge, this theoretical work has not been ap- 
plied to any real case yet. A basic and simple idea to 
take vibrations into account is to calculate the absorp- 
tion cross-section for a configuration where the absorbing 
atom is shifted from its equilibrium position^i^ For in- 
stance, at the B K-edge of ABO3 perovskite structures, 
the approach used by Vedrinskii et al^ to treat the tem- 
perature effects is based on the assumption that the area 
of the pre-edge peak is proportional to the mean square 
displacements of the B atom along the 0-B-O chain. In 
such a way, the symmetry-breaking of the absorbing atom 



site generates or increases local p-d and p-s hybridiza- 
tion. Brouder and co-workersii5i developed a different 
approach assuming that vibrational energies are small 
with respect to the instrumental resolution, and using 
the crude Born-Oppenhcimer approximation so that only 
core-hole motion remains. The resulting expression for 
the absorption cross section shows that, at the K edge, 
vibrations enable electric-dipole transitions to 3s and 3d 
final states, that are not due to local hybridization with 
the p states. In the case of organic molecules, more ad- 
vanced theories are employed to reproduce the vibronic 
structure observed in XANES spectra. Beyond the Born- 
Oppenheimer approximation, a vibronic coupling the- 
ory was developed and successfully applied to the C- 
Is absorption spectra of ethylene (C2H4)^^ and ethyne 
(C2H2)^^ and isotopomers. Besides, many approaches 
calculate Franck-Condon factors based on a vibrational 
eigenmode analysis in the ground and excited states. In 
particular, the linear coupling model was used to com- 
pute the Franck-Condon factors for NEXAFS spectra 
at the O K-edge in formaldehyde^^ at the C iiT-edge 
in formaldehyde^^ naphtalene^, gaseous cyclopropane^ 
and acetic acid^^ at the N K-edge in gaseous pyridine^^ 
acetonitrile^ and acrylonitrilcj^ In Ref. |9| a direct calcula- 
tion of the Franck-Condon amplitudes is also performed. 
However, the Franck-Condon approximation ignores the 
impact of nuclear motion on the electronic transition 
amplitude. To first order, this impact is referred to as 
the Herzberg- Teller effect.^'* To go beyond the Franck- 
Condon approximation, the nuclear degrees of freedom 
of molecules was modeled using various molecular dy- 
namics (MD) techniques. In doing so, atomic configura- 
tions are generated at finite temperature and then used 
as input in cross-section calculations. Finally, the indi- 
vidual spectra associated to each sampled atomic con- 
figuration are averaged over. Using this methodology, 
vibrations were included in XANES calculations at the 
N K-edge of several prototype molecules: upon classical 
MD samplingji^^— and upon path-integral molecular dy- 
namics (PIMD).— i^ PIMD sampling enables to consider 
the quantum nature of nuclear motion, which has been 
found to be of noticeable importance to accurately sim- 
ulate XANES of N2iS£, and s-triazine and glycine^ The 
sensitivity of photoabsorption spectroscopy to a quantum 
treatment of nuclear motion was also highlighted in the 
optical ranged and in the UV range»22. The MD method- 
ology is not restricted to organic molecules; simulations 
of the density and temperature dependence of XANES 
in warm dense aluminum plasmas were achieved using 
ab initio2^2^ and classical^^ MD and recently they have 
been compared to experiments i^iiZ^ 

In the present study, DFT calculations based on plane- 
wave formalism are performed in order to understand 
the spectral modifications induced by temperature that 
are observed experimentally. Reference theoretical Al 
K XANES spectra of corundum and beryl are first ob- 
tained considering the atoms at their equilibrium posi- 
tions. Then three existing methods are used to account 



for vibrations: (i) a calculation considering thermal dis- 
placements of the absorbing atom around its equilibrium 
position, (a) the method developed by Brouder et al.^ 
(in) a convolution of the spectra obtained for the atoms 
at the equilibrium positions with an approximate phonon 
spectral function. 

The paper is organized as follows. Section HI] is dedi- 
cated to the methods. First the temperature-dependent 
XANES experiments at the Al K-edge of corundum and 
beryl are described, with a specific attention paid to 
the self-absorption correction procedure. Second, the 
three theoretical DFT methods used to take the vibra- 
tions into account are detailed. Section HTT] is devoted to 
the description of the experimental spectra recorded for 
temperatures ranging from BOOK to 930K. In section ITVl 
the experimental results are analyzed with the help of 
the DFT calculations, and discussed in the context of 
the temperature-dependent XANES spectra already re- 
ported in the literature. The conclusion of this work is 
given in Section [Vj 



II. METHODS 
A. Experimental setup 

Two single-crystalline samples containing aluminum in 
6-fold coordination with oxygen are studied. The first is 
a cylindrical synthetic transparent ruby (a— AI2O3 with 
18 ppm of Cr and 10 ppm of Ti) of 4.0 mm diameter 
and 1.0 mm thickness, which will thereafter be denoted 
as corundum. The second is a parallelepipedic section of 
a natural green emerald from Colombia (Be3Al2Si60i8 
with 238 ppm of Fe, 78 ppm of Cr, 59 ppm of Sc and 
36 ppm of Co) of 10.0 mm length, 3.6 mm width and 
1.2 mm thickness, which will be denoted as beryl. The 
samples were analyzed using the CAMPARIS electronic 
microprobe at Universite Pierre et Marie Curie - Paris 6, 
France. 

X-ray absorption experiments were performed at 
LUCIA beamline in the French synchrotron facility 
SOLEIL.^^ The synchrotron was operating in the top-up 
mode at 2.75 GeV with a current of 400 mA. The beam 
spot size was set to 1 x 2 mm^ and the energy range cho- 
sen to include the Al K-edge (1559.6 eV in bulk)^^ was 
set from 1550 eV to 1700 eV, hence measuring pre-edge, 
XANES and the beginning of EXAFS. Energy selection 
was performed through the combination of an HU52 "Ap- 
ple II" type undulator and a double KTP (Oil) crystal 
monochromator. The pressure in the experimental cham- 
ber was 10^^ mbar. 

The samples were held between a parallelepipedic 
boron nitride furnace and a punched molybdenum 
lamella fixed to the furnace, allowing temperatures rang- 
ing from 300K to 950K. This holder was fixed on a mobile 
stage, allowing translations along a cartesian xyz coor- 
dinate system plus a 2-axis rotation, where the cc-axis 
corresponds to the beam direction and the y-axis to the 



(horizontal) linear polarization e of the beam. In order 
to minimize self-absorption effects, the sample surfaces 
should be orthogonal to the beam axis. For fluorescence 
detection, the sample holder has been slightly rotated by 
15° around z-axis. 

Corundum and beryl point groups are 3m and —mm 
respectively] ^^1^^ This implies that both materials are 
dichroic in the electric dipole approximation. ^° The 
single- crystals were set on the sample holder in order 
to measure successively a\\ and a± spectra, correspond- 
ing respectively to e parallel and perpendicular to the 
high symmetry axis, i.e. the three-fold symmetry axis 
for corundum and the six-fold symmetry axis for beryl. 

Total fluorescence yield was measured by a four ele- 
ment silicon drift diode (SDD) detector with a total ac- 
tive area of 40 mm^ protected from infrared and visi- 
ble radiations by a thin beryllium window. In order to 
maximize the signal/noise ratio, each point was obtained 
after a six second acquisition time and five consecutive 
spectra were measured for each configuration (a config- 
uration consists of a sample, an orientation and a tem- 
perature). This system produces five intensity outputs, 
one for each of the four SDD (/i, I2, I3, h) and one 
for the incident beam (/q); measured before the sample. 
For each spectrum, each Ij was then divided by Iq, fol- 
lowed with a normalization on the 1590-1700 eV energy 
range (spectra are divided by their mean values on this 
region) giving (/,//o)norm- Then, (/,//o)av was calcu- 
lated for each configuration by averaging five consecutive 
acquisitions of (/j//o)norm- Since self-absorption effects 
are angle dependent and since the four SDD were not at 
the exact same spatial position, the four (/j//o)av were 
not directly comparable. Self- absorption correction was 
performed using the formulai ^^'^^ 
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j, Nj{uj) — (J ij)}^ \ is the measured intensity 
on silicon drift diode j, normalized at a fixed pulsa- 
tion cjQ (chosen after the main edge at ?iajo=1653 eV), 
6i=75° is the angle between the beam axis and the sam- 
ple surface, 9fj is the angle between the sample sur- 
face and the outgoing fluorescence beam (which depends 
on each detector j position). Coefficient f3 is equal to 
^*°\cj")'° where Wfluo is the Kai-fluorescence pulsation 
(?icjfluo=1486.7 eV) and /^tot is the total absorption coef- 
ficient, including edge absorption /i and background con- 
tribution /ibg of other atomic species in the sample. Fi- 
nally 7 is equal to ^^ (wo ) . Values for /3 and 7 are calcu- 
lated based on crystal stoichiometry, atomic masses, unit 
cell geometry and using experimental tabulated atomic 
values,^ For example, Mtot(w)=X]at /^atAil?'''°(w)/Mat 
with Pe,t=Na,tMs^tJ^A/Vu.c. whcrc the sum is over all the 
atoms in the unit cell and for each atom /i*t'''° is the 
tabulated value in Barns/atom, Mat is the molar mass 



in g.mol"^, TVat is the number of atoms 'at' in the unit 
cell, Vu.c. is the unit cell volume in cm^ and AFa is the 
Avogadro constant. Such a correction gives comparable 
values of aj for each detector, each one being normalized 
at the point of energy hajQ. By using such a procedure 
to perform self-absorption correction, averaging over the 
four detectors is possible and gives the spectra a\\ and 
a± shown in section IIIII No smoothing were applied to 
the spectra. 



B. Theoretical methods 

In order to explain the experimental spectral features 
observed when temperature is increased, first-principles 
calculations were performed in the density functional 
theory framework using the Quantum-Espresso suite of 
codesi^ The Is core-hole effects were taken into account 
within a 2x2x2 trigonal supercell for corundum, con- 
taining 80 aXoTasr^ and a 2x 1 x 1 hexagonal supercell for 
beryl, containing 58 atoms. ^^ Troullier-Martins^^ norm- 
conserving pseudopotential were used. The pseudopoten- 
tials of Al and Si were generated using the 3s, 3p and 3d 
orbitals as valence states, with 2.00 Bohr cutoff radii and 
considering the d states as local. The pseudopotential of 
the Al absorbing atom was generated with only one Is 
electron in the Al electronic configuration. The oxygen 
and beryllium pseudopotentials were built using the 2s 
and 2p as valence states (local part p) with cutoff radii of 
1.46 Bohr for O and of 2.00 Bohr for Be. Self-consistent 
charge densities were calculated in the generalized gra- 
dient approximation of ReflSg. using the PWscf code, at 
the r point, with plane- wave energy cut-off of 80 Ry. 
Thereafter, local and partial density of states were cal- 
culated using Lowdin projections on a 4x4x4 fc-point 
grid with a Gaussian broadening parameter of 0.3 eV. 
XANES cross section calculations were performed using 
the XSpectra code .^^i^^ The all-electron wave-function 
of the Al absorbing atom was reconstructed using the 
Projector Augmented Wave method,*^ with augmenta- 
tion region radii of 2.0 Bohr. The spectra were com- 
puted on 4x4x4 fc-point grid with a broadening param- 
eter 7=0.6 eV. 

Three methods were used to take into account the ef- 
fects of vibrations without having to calculate phonon 
modes. First, adapting the idea of RefUSl on vibra- 
tions, only the absorbing atom was moved in the crystal 
structure with respect to its equilibrium position. This 
method will be referred to as 'method 1' in the follow- 
ing section. As vibrations are concerned, spectra for sev- 
eral absorbing atom displacements R were calculated and 
averaged with a weighting core displacement distribution 
function p(R) . Within the harmonic approximation, this 
distribution function is linked to the thermal parameter 
matrix [Uij] through p^R) = exp(— R • ^^ • R). The 
components of the rank-2 tensor U are given by X-ray or 
neutron scattering measurements. In this approximation 
each spectrum is calculated with a static off-center ab- 






FIG. 1. A sketch of the 21-point cubic grid used in Eq. (2] 



sorbing atom, creating hybridization between Al valence 
orbitals and neighboring atoms orbitals, thus modifying 
valence states of the crystal. Vibration modeling comes 
from averaging spectra over a large number of absorbing 
atom displacements. 

The absorbing atom motion was discretized on a cubic 
tridimensional grid consisting of 361 points. Equation [2] 
gives the cubic integral approximation on a grid of 21 
points (see Fig. [T]):^^ 
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(2) 
where a is the length of the cube edge, 
/(a;,R)=CTe(aj,R) x /9(R) is the weighted spectrum 
calculated with absorbing atom displaced of R, 
frn = /(w,0) is the value at the center of the cube 
(red point), ^ fr is the sum of the values of / at the 4 
points halfway between the center of the cube and the 
center of each face (pink points), ^ // is the sum of the 
values of / at the center of each face (green points) and 
X] fd is the sum of the values of / at the vertices of the 
cube (blue points). The 361 point grid is an augmented 
3x3x3 version of the 21 point one. This discretization 
would lead to a 27x21 point grid, which reduces to 
only 361 points after avoiding double counting of similar 
points. The value of the integration cube edge a is 
given by 2^JUeig^ where f/cig are the eigen values of 
the U matrix. For the 'method V calculations, a was 
chosen equal to 0.3 Bohr radius, in agreement with the 
thermal matrix parameters of corundum and beryl given 

inRefs.|7§and|9l!- 

The second method will be referred to as 'method 
2'^i56 gjj^(3g vibrational energies are small with respect 
to the instrumental resolution, the closure approxima- 
tion is made. In this approximation, the final vibrational 
states of the electrons+nuclei system are summed over. 
Hence atomic motions are assumed not to change the 
final electronic states that are calculated at the atomic 
equilibrium positions. Consequently, in 'method 2', only 
the position of the absorbing atom in the initial state is 
modified by vibrations. In practice, for a K-edge, the Is 
wave-function of the initial state is no more centered on 



the absorbing atom equilibrium position. This shifted Is 
wave-function is expanded over spherical harmonics cen- 
tered on the equilibrium position. The H.—1 component 
allows the transition to 3s states. Contrary to 'method 
1', this method does not create hybridization in the fi- 
nal state. Calculating absorption cross section in this 
method is based on Eq. (9) from Ref. [l]. The Is wave- 
function movements are calculated using the same core 
displacement distribution used in 'method 1', integrating 
over the same 361 point grid, but using a shorter value 
for a (a=0.1 Bohr radius as in Ref. [1]). 

The third method, denoted as 'method 3', is based 
on Fujikawa's work/^ Within the Franck-Condon ap- 
proximation, the effect of vibrations on x-ray absorp- 
tion spectroscopy can be represented as the convolution 
of the phonon spectral function with the x-ray absorp- 
tion spectrum at equilibrium position. Here the phonon 
spectral function is approximated as a Gaussian distribu- 
tion. Along with the two precedent methods, convoluted 
spectra are calculated. The convolution is performed 
between the first-principles equilibrium spectrum for a 
configuration and a Gaussian. The Gaussian equation is 
l/(-\/27rcr) exp(— x^/(2cr)^), with parameter a being cho- 
sen for each configuration in order to reproduce the main 
edge intensity of method 1 spectra. 

Calculated spectra presented in section HVl are all nor- 
malized on the high-energy region, as was done for ex- 
perimental data. 



III. EXPERIMENTAL RESULTS 

Figure [D shows the self-absorption corrected XANES 
spectra of corundum and beryl recorded at different tem- 
peratures for both orientations. Spectra originating at 
2/=— 0.5 correspond to the difference (t(T)— (7(300K). We 
observe here for the first time the slow evolution of the 
Al K pre-edge feature with temperature. In addition to 
the expected increase in the peak intensityjii^ a shift to- 
ward lower energy of the pre-edge structure is observed 
when temperature increases. Furthermore, this effect ap- 
pears for each sample and for each orientation and seems 
to be general, as well as the decreasing intensities and 
broadening of main peaks with temperature. 

Corundum (T|| and a±_ spectra are quite similar up to 
1575 eV: they both present a well-resolved pre-edge peak 
P, the main peak A followed by the double feature Bi and 
B2. The 300K polarized spectra are in agreement with 
those recorded in total electron yield shown in Ref.lO^. As 
mentioned before, the intensity of peak P increases with 
temperature while its position is shifted towards lower 
energy. This is confirmed by the peak p in the difference 
spectra. Peak A also varies with temperature: it broad- 
ens, its intensity decreases and its position is shifted to 
lower energy (around 0.2 eV and 0.3 eV for a\\ and crj_ re- 
spectively). The decrease of peak A intensity corresponds 
to the difference spectra peak a and the shift to lower en- 
ergy to a ' . Peak a ' is also related to the broadening of 
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FIG. 2. Experimental X-ray absorption polarized spectra of corundum (left) and beryl (right) at the Al A'-edge for different 
temperatures along with difference of each spectrum with respect to the 300K reference. The top and bottom panels display 
the CTii and (tj_ spectra, respectively. 



peak A with temperature. Peaks Bi and B2 intensities also 
diminish and they only form one broad peak at 930K in 

cr|| and a±. 

Beryl a\\ and a± spectra present much different shapes, 
unlike the case of corundum. In the energy range 1560- 
1575 eV, the (T|| spectrum exhibits the pre-edge peak P 
and three main features A, B and C, feature B being around 
twice more intense than A and C at 300K. On the con- 



trary, the CTj_ spectrum exhibits two main features B and 
C, preceded by a small shoulder A and the pre-edge P. 
Peak C is more intense than peak B. Peak A and B of a± 
are located at 0.7 eV and 0.3 eV lower than their corre- 
sponding peaks in crii . The pre-edge peak P is less intense 
in (7j_ than in a\\ at 300K and its intensity increases more 
in cr|| than in a±. Indeed the difference peak p appears 
better resolved and more intense in a\i than in a±. Peak 



A of (T|| has a constant intensity with temperature but cre- 
ates the difference peak a ' due to broadening. The weak 
peak A of <t^ decreases in intensity with temperature. 
The B peak intensity also decreases with temperature, 
as well as a slight shift of 0.1 eV toward lower energy, 
that creates peaks b and b' in the difference spectra, 
more pronounced in crii than in a±. Peak C decreases 
and broadens with temperature, leading to feature c in 
the difference spectra. 
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apparently a key feature of the temperature evolution 
of the spectra. The features at high energies are only 
slightly modified. Sum rules relates XNLD to the electric 
quadrupole distribution of the empty DOS with p sym- 
metry on the aluminium site.^'^ The sum rules were first 
derived in the absence of electron-phonon coupling and 
it is not yet clear whether the sum rules still hold when 
temperature is present. A theoretical analysis of the tem- 
perature dependence of XNLD sum rules is beyond the 
scope of the present paper but would eventually give in- 
formation on the modification of the electronic structure 
with the temperature. 

In order to more quantitatively describe the behavior 
of the pre-edge peak, extraction was performed using Fi- 
tyk software,—- which fits the main edge and pre-edge by 
the sum of an arctangent and a Gaussian. Figure [5] re- 
ports the fitted Gaussian intensities and center energies 
for each sample and each orientation. These quantitative 
values are in good agreement with the foregoing obser- 
vations, a global increase of intensity with temperature 
is noticeable. While both orientations seem to progress 
similarly in corundum, it is not the case in beryl where 
a± intensities grow less rapidly than in cry . Concerning 
peak position, both compounds and orientations present 
a similar shift toward lower energies. 



IV. INTERPRETATION AND DISCUSSION 

In this section, theoretical results are presented and 
discussed in order to try to understand the effects of tem- 
perature observed in the XANES spectra shown in Sec- 
tion [nil First, a density of states (DOS) analysis is per- 
formed for both compounds showing the involvement of 
the 3s empty states of the Al absorbing atom in the pro- 
cess of pre-edge peak creation. Then, theoretical spectra 
obtained by using the three different methods described 
in section fll Bl are shown and compared together. 



FIG. 3. Temperature dependence of the X-ray natural lin- 
ear dichroism (XNLD) measured at the Al ic'-edge in corun- 
dum (top panel) and beryl (bottom panel). The cross section 
o"XNLD is equal to the difference a^^-a±. Note that the j/-axis 
range is twice larger for beryl than for corundum. 



The x-ray natural linear dichroisms (XNLD) for the 
various experimental temperatures are plotted in Fig. [S] 
The XNLD signals of the two samples are large and their 
variations with temperature do not change their general 
shape. XNLD on the pre-edge is quite small so that tem- 
perature dependence can hardly be seen. One mainly 
observes that some sharp XNLD maxima are reduced 
when the temperature increases: Al and A2 peaks for 
corundum and B peak for beryl. One also notices that 
the rising edge and the maximum of the first XNLD fea- 
ture is shifted to lower energies. This effect is observed 
for both corundum and beryl and can also be detected 
on the parallel and perpendicular cross-sections. It is 



A. DOS calculations: the nature of the pre-edge 
peak 

The pre-edge peak is not reproduced by purely elec- 
tric dipole Is— J-p transitions, as shown in the top panels 
of Fig. \5\ Hence some transition toward non-p states 
must be at stake. The bottom panels of Fig. [5] present 
partial and local density of states (DOS) of corundum 
and beryl. The s, p and d partial empty DOS are plot- 
ted for the absorbing aluminum atom (with a core hole) , 
only the s and p for the first six oxygen neighbors, and 
again the s, p and d for the aluminum next neighbor. 
The similarity between the absorbing Al 3p local DOS 
and the XANES spectra shows that the later is a good 
probe of the former. In the pre-edge region, the role 
of the absorbing Al atom is clear as the pre-edge po- 
sition coincides with its 3s projected DOS. In the case 
of beryl, a certain proportion of absorbing Al 3d states 
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FIG. 4. Pre-edge intensity and energy position for both orientations and for both samples as a function of temperature. The 
y axis units of the upper panels are consistent with those of Fig. [2] 



is also present in the pre-edge region. However, the 3d 
states do not contribute to the absorption coefficient via 
electric quadrupole (Is— ^3d) transitions. Indeed the cal- 
culated electric quadrupole transitions are negligible in 
the whole XANES region, since they are found to be 
lower than 0.25 10~^ in the normalized units used in the 
top panels of Fig. [5] The 3s states of the absorber are 
hybridized with both s and p empty states of the oxy- 
gen neighbors while the aluminum next neighbors do not 
contribute to the pre-edge since their DOS are at higher 
energies. The presence of the Is core hole on the absorb- 
ing atom shifts the first empty states to the bottom of 
the conduction band and leads to the pre-edge and first 
peaks of XANES. Non-excited aluminum atoms do not 
contribute to the pre-edge feature of the spectrum. 



B. Theoretical XANES spectra: vibration 
modeling 

The DOS analysis indicates that the modeling of the 
experimental Al i^-edge XANES spectra has to allow 
transitions toward the 3s states of the absorbing Al. Both 
method 1 and method 2 presented in section IIIBI are 
able to achieve this goal. Figure |6] displays the theo- 
retical spectra obtained by the three methods described 
in section Hi Bl together with the spectra calculated with 
the atoms at their equilibrium positions (labeled 'equil.'). 
Difference spectra with respect to the 'equil.' reference 
are also plotted. It should be noted that the 'equil.' spec- 
tra do not strictly correspond to the OK case since they 
do not take the zero point motion of nuclei into account. 
Experimental spectra measured at 20K on the same sam- 
ples (not shown here) still exhibit the pre-edge peak, even 



without active phonon modes, and the pre-edge intensity 
at 20K is quite similar to that at 300K. At room temper- 
ature nearly no phonon modes are activCj-^'-^ thus the 
pre-edge feature observable in the 20K-300K range seems 
to be essentially due to the quantum zero point motion 
effect. 

Corundum aj_ and <t\\ 'equil.' spectra reproduce well 
the overall shape of the experimental spectra, except in 
the pre-edge-region. Switching on method 1 mostly af- 
fects the main edge A intensity and shifts its position by 
0.2 eV with respect to 'equil.' spectra. Peaks a and 
a' of the difference spectra are the signature of these 
effects. Peak Bi intensity slightly decreases for both ori- 
entations while peak B2 is nearly unchanged. Small dif- 
ference peaks bl and b2 confirm these points. A weak 
pre-edge peak P seems to appear for both orientations. 
Although method 1 is not able to grow a well-resolved 
pre-edge peak, it provides difference spectra that are very 
similar to the temperature-dependent experimental dif- 
ference spectra of Fig. IIIII (left panels). On the contrary, 
method 2 gives rise to a more important pre-edge peak, 
which is in better agreement with experiment. The differ- 
ence peak p shows a really distinct Gaussian-like shape. 
The rest of the spectrum is nearly unchanged and this 
is confirmed by the flat difference spectra. Method 3 
spectra do not show a pre-edge but have a visible im- 
pact on main edge peak A, by broadening it, lowering 
its intensity, and shifting its energy to higher energies 
by 0.1 eV, i.e. in the opposite direction as compared to 
temperature-dependent experiments. These effects are 
visible on difference peaks a and a ' . 

Beryl an and cr± 'equil.' spectra also reproduce the ex- 
perimental spectra quite well, except for some of the peak 
intensities. While peaks A and C are of similar amplitude 
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FIG. 5. Top: Calculated and experimental isotropic (2a"x+o"||)/3 spectra at the Al K-edge of corundum (left) and beryl(right). 
Bottom: partial and local density of states of corundum (left) and beryl (right) calculated for the supercell including a Is core 
hole. The role of the Al 3s of the absorbing atom is prominent in the pre-edge region. 



in CT|| experimental spectrum, the intensity of peak A is 
twice the one of peak C in the 'equil.' spectrum. On 
the opposite, peaks B and C have similar intensity in a±_ 
'equil.' spectrum while peak B is 1.5 times less intense 
than peak C in the 300K experimental spectrum. By 
using method 1 for (T|| , no pre-edge peak seems to arise, 
peak A is decreased and shifted to lower energy by 0.2 eV, 
peak B is also decreased and peak C remains stable. These 
effects are also visible through the shape of the difference 
spectrum. Thus vibrations treated within method 1 are 
not able to well reproduce the temperature-dependence 
trend observed in the top-right panel of Fig. [21 In the 
case of a±, method 1 leads to theoretical spectrum in 
better agreement with the x-ray absorption temperature- 
dependence seen in the bottom-right panel of Fig. [2] In- 
deed, method 1 enables a well resolved pre-edge peak P 
to rise (clearly confirmed by difference peak p), and a 
very weak peak A to appear. It also moves B to lower en- 
ergy (0.2 eV) and decreases the intensity of peaks B and 
C. Here again, method 2 has no influence on peak posi- 
tion, it nearly does not change the amplitude of peaks 
A, B, C in cr|| and B, C in (t±. However it creates the 
pre-edge P at the right energy, but with too high an in- 
tensity. The value of the U parameter used in method 2 



might be overestimated for beryl and yielded a too pro- 
nounced pre-edge peak. Method 3 decreases the ampli- 
tude of main peaks. It also shifts peak A of cr|| and B of 
CTx towards higher energies that, as in corundum, is the 
opposite of what is seen experimentally. A slight broad- 
ening appears and is more visible in a±. 

To summarize, method 1 leads to difference spectra 
that reproduce quite well the temperature-dependent 
trend observed experimentally in most cases. Neverthe- 
less the impact of vibrations as treated by method 1 is 
not sufficient to systematically reproduce the experimen- 
tal pre-edge feature. The p-s hybridization induced by 
this method is not satisfactory enough. In method 1, 
the vibration modeling results in integration, over a cu- 
bic volume around the absorbing atom, of cross sections 
calculated for absorbing Al isotropically displaced inside 
this volume. In other words, this method equally consid- 
ers the contributions of Al motion along all directions. 
Some Al displacements would be favored by taking into 
account the phonon modes of the material, what could 
substantially improve the agreement between theoretical 
and experimental spectra at least in the pre-edge region. 
Method 2 spectra all exhibit a well-defined pre-edge peak, 
corresponding to dipole ls-3s transitions, which are al- 
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FIG. 6. Calculated a\\ and a± spectra at the Al K edge of corundum (left) and beryl (right) along with difference of each 
spectrum with respect to the equil. reference. The 'equil.' spectra correspond to DFT calculation performed with the atoms 
at the equilibrium positions. The spectra labelled 'meth. 1', 'meth. 2' and 'meth. 3' refer to the three methods of vibration 
modeling described in section [HI The upper-right inserts use the same legend code and present a zoom on the main peak(s). 



lowed by the displacement of the Is wave-functioniii^ 
Nevertheless, method 2 has nearly no impact on the peak 
intensity above the edge, and is unable to produce any 
shift of the features. This last result was expected be- 
cause of the use of the crude Born-Oppenheimer (BO) 
in method 2. Within the crude BO, the transition en- 
ergies are evaluated for the atoms at their equilibrium 



positions, hence the spectral features remain at the same 
energy positions as in 'equil.' spectra. Method 3 spec- 
tra, as expected, do not exhibit any pre-edge peak, and 
present main-edge peaks less intense and shifted in en- 
ergy. However the energy shift goes the wrong way, i.e. 
to higher energies. Therefore, a simple convolution of 
the 'equil.' spectra seems to be inappropriate to account 
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for the temperature-dependence of the XANES spectra 
observed experimentally. 

Thermal expansion effects have been evaluated in the 
case of corundum by performing a XANES calculation 
for the structural parameters refined at 2170K (153K be- 
low the melting point).™ A contraction of the highest 
energy XANES features has been noticed. Such an effect 
was expected according to the predictions of the Natoli's 
rulei^ However, it is contradictory with what has been 
experimentally observed when temperature is increased 
(Fig. [2|). Therefore, the thermal expansion is unable to 
explain the temperature-dependence of the Al K-edge 
XANES spectra in corundum and beryl. 

The calculations carried out in this study clearly show 
the crucial role of vibrations in the pre-edge region and 
partially explain the temperature-dependence observed 
in the whole XANES region. However, none of the three 
methods used here is fully satisfactory. A great im- 
provement would be to use PIMD as Schwartz et al. 
did in the case of two isolated organic moleculesi^ But 
the computational cost of PIMD simulations in solids, 
such as corundum or beryl, would be a serious limit- 
ing factor. A more reasonable way to account for the 
thermal fiuctuation of XANES at the Al K-edge would 
consist of the generation of atomic configurations of the 
whole crystal from the dynamical matrix of the system. 
The temperature-dependent XANES theoretical spec- 
trum would then result in an average of cross-sections 
calculated for a large number of configurations. An anal- 
ogous methodology has been successfully employed to ac- 
count for the temperature-dependence of nuclear mag- 
netic resonance chemical shift in MgO<^ 



V. CONCLUSION 

In this study, temperature-dependent Al K-edge 
XANES spectra of corundum and beryl have been pre- 
sented for temperature ranging from 300K to 930K for 
the first time. The XANES spectra were measured on 
single- crystals with the polarization vector of the x-ray 
beam parallel and perpendicular to the high symmetry 
axis of the minerals, allowing the investigation of the tem- 
perature dependence of the CT|| and <7± components. This 
series of experiments shows that the pre-edge peak is very 
sensitive to thermal fluctuations: the pre-edge peak in- 



tensity grows and its position shifts to lower energy as 
temperature increases. These variations do not depend 
much on the x-ray polarization. Thermal fluctuations 
are also visible above the pre-edge region, through an in- 
tensity decrease of the main features and, in the case of 
corundum, through a slight shift to lower energy of the 
flrst main peak. 

First-principles DFT calculations have conflrmed the 
Al 3s nature of the pre-edge peak in both minerals, and 
have shown that the introduction of vibrations within the 
XANES calculation gives rise to a pre-edge peak. Our 
experimental and theoretical results on the Al K pre- 
edge in corundum and beryl suggest that vibrations (the 
zero point motion and also temperature) could be able 
to induce a pre-edge feature at the Al K-edge in other 
minerals. This conclusion might bring a reinterpretation 
of the appearance of a pre-edge peak observed at the Al 
K-edge in zeolites when temperature is increased, which 
was interpreted as the signature of three-fold coordinated 
Al.i^iiiii^ 

We have used three existing theoretical methods to 
take nuclear motion into account. These methods have 
provided useful information to understand the origin of 
the spectral modiflcations observed in the temperature- 
dependent experimental XANES spectra. However, they 
were found to be inappropriate to fully account for 
the thermal fluctuation of XANES. This knowledge will 
be particularly useful for upcoming works, giving in- 
sight of what has to be considered in the modeling of 
temperature-dependent XANES spectra. 
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